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Abstract 



The quantum effective Yang-Mills condensate (YMC) dark energy 
model has some distinguished features that it naturally solves the co- 
incidence problem and, at the same time, is able to give an equation 
of state w crossing —1. In this work we further employ the Statefinder 
pair (r, s) introduced by Sahni et al to diagnose the YMC model for 
three cases: the non-coupling, the YMC decaying into matter only, 
and the YMC decaying into both matter and radiation. The tra- 
jectories (r, s) and (r, g), and the evolutions r{z), s{z) are explicitly 
presented. It is found that, the YMC model in all three cases has 
r ~ 1 for z < 10 and s ~ for z < 5 with only small deviations 
~ 0.02, quite close to the cosmological constant model (LCDM), but 
is obviously differentiated from other dark energy models, such as 
quiesence, kinessence etc. 

keywords: Yang-Mills condensate, dark energy, Statefinder 
PACS numbers: 98.80.-k, 98.80.Es, 95.36.+X, 04.62.+V 



*mltong@mail.ustc.edu.cn 



1 Introduction 



The currently accelerating expansion of the universe has been indicated by the observation of 
Type la supernovae [TJ, and is consistent with the data from cosmic microwave background 
(CMB) [21 [3] and from the cosmic large scale structure [3]. The acceleration of the expanding 
universe is attributed to the dominant dark energy Qj>^ ~ 0.73, coexisting with the matter 
Qm — 0.27. The simplest model for dark energy is the cosmological constant [5J, with an 
equation of state (EOS) w = —1 as the universe evolves. However, two questions arise from 
this scenarios, namely, the fine-tuning problem and the cosmic coincidence problem. While 
the former exists for almost all the dark energy models, the latter has been addressed in the 
class of dynamical dark energy models, which take some dynamically-evolving field as the 
candidate for the dark energy. Among them are the quintessence [6], phantom [7], k-essence 
[8], quintom [9], tachyonic [lOj, holographic dark energy model [IT], and interacting dark 
energy model [12]. Besides, there is another interesting type of model built on the quantum 
effective gravity [13] , which includes the quantum corrections of gravitational field to Einstein 
equations. Different from these models, the Yang-Mills condensate dark energy model is based 
on a vector- type of the quantum effective Yang-Mills field [HJ [151 [IS [13 [Hj- From field- 
theoretical point of view, the model has the following interesting properties: the gauge fields 
are indispensable to the Standard Model of particle physics, and the effective Lagrangian of 
YMC is determined from the standard field-theoretical calculations for each order of loops of 
quantum corrections, and thereby gives the correct trace-anomaly, and there is room for change 
its form by hand. Moreover, it is found that, for quite generic initial conditions, the YMC dark 
energy model always has the desired tracking behavior that naturally solves the coincidence 
problem. This has been accomplished for the cases of 1-loop [I6lll7|, 2- loop |18j, and 3-loop 
|19j quantum corrections, either with coupling or without coupling to matter and/or radiation. 
When coupling with matter, or with both matter and radiation, the YMC has an EOS crossing 
— 1 smoothly and taking w ~ —1.05, as indicated by the recent preliminary observational data 
of the Supernova Legacy Survey [201 I^T]. 

In order to differentiate these various dark energy models, Sahni et al [22] introduce a 
new geometrical diagnostic pair (r, s), called Statefinder, which involves the third order time- 
derivative of scale factor. The pair is related to the EOS of dark energy and its time derivative. 
From the observational side, the values of (r, s) can be extracted from data coming form 
SNAP type experiments [23]. The Statefinder diagnosis has been applied to several dark 
energy models [211 [23 [23 [2Z]. In particular, the spatially flat LCDM has a fixed point 
(r, s) = (1, 0). For the 1-loop YMC dark energy model, Ref.[28] studies the non-coupling case 
with the radiation contribution being neglected. In this paper for a complete treatment of the 
Statefinder diagnosis, we work with the 2-loop coupling YMC model and include the radiation 
component. Thereby, there arise considerable modifications to the 1-loop non-coupling case. 



2 2-loop YMC dark energy model 



We consider a spatially flat (A; = 0) Robertson- Walker (RW) universe, whose expansion is 
determined by the Friedmann equations 
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where H = ^, the pressure p(t) = Py + Pr, the energy density p(t) = Py + Pm + Pr, with the 
subscripts 'y', 'm' and 'r' refer to the YMC dark energy, the matter (including both baryons 
and dark matter), and the radiation, respectively. Up to the 2-loop quantum corrections, the 
energy density py and the pressure Py of the YMC are given by [HI [19] 
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where y = In \F/k'^\, F = E"^ — B^, k is the renormalization scale with dimension of squared 
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representing the 2-loop contribution, and the dimensionless constant 5 is a parameter 
representing higher order corrections. For simplicity, we take the gauge group to be SU{2) 
and only consider the 'electric' condensate, i.e., F = E"^ . The EOS for the YMC is 
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When one sets r/ = in the above expressions of Eqs.Q, (jlj), and ([5]), the 1-loop model is 
recovered [IHl [lEl |l7j . At high energies y oc, Py is positive, and the EOS of YMC approaches 
to that of a radiation w 1/3, as is expected for an effective quantum field theory. 
The dynamical evolutions of the three components of the universe are 



Py + 3-{py+ Py) = -Tpy " F'^J/, 
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pr + 3-{pr+Pr) = ^'py, 



where F and F' are the decay rate of YMC into matter and radiation, respectively. Taking py in 
Eq.(l3]) to be the dark energy with Py{to) ~ 0.73pc, the scale k is fixed by k^/^ ~ l.QhH"^ x 10~^ 
eV, where /iq ~ 0.72 is the current value of the Hubble parameter. For concreteness, 5 = 3 
and 7] = 0.84 are taken. 

The initial matter and radiation densities at z ~ 3454 [3] are taken to be = -Rj ~ 1 x 10^°, 
where x = pm/ \bK^, and R = pr/ The initial YMC can be chosen to be in a broad range 
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Figure 1: The evolution of tlie energy densities. In the non-coupling case pm oa a{t)~^ and pr oc a(i)~'* 
throughout. In the case of YMC decaying into matter, p,„ deviates from oc a{t)~^ around z ~ and levels 
off. In the case of YMC decaying into both matter and radiation, both p„i{t) and Pr{t) level off around z ~ 0. 
Note that two different initial condition: flyi = 10^^ and 10^^", yield the same py, p,„, and pr at z ^ 0. 

Ui = (1, 15), corresponding to Qyi ~ (10~^°, 10~^). For the case of YMC decaying into matter, 
we take T / Hq = 0.2 with Hq = approximately the Hubble constant. For the case 

of YMC decaying into both matter and radiation, we take T/Ho = 0.2 and T' /Ho = 1.2 x 10~^. 
Fig{T] shows the evolution of the energy densities for various components in all three cases. It 
is seen that, for the initial values of Qyi ranging eight orders of magnitude, the present status 
Qy ~ 0.73, Qm ~ 0.27, and Qr ~ 10~^, are achieved. In this sense, the coincidence problem 
of dark energy is naturally solved in the YMC model [151 [13 [IB]- K is also found that in the 
coupling cases w can be cross —1 smoothly [[18] . 

3 Statefinder 

The Statefinder pair (r, s) are defined as the following [22] 
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where the deceleration parameter 

q = — '^ = hi + 3nyw + nr), (10) 

with Qy = Py/p and Qr = Pr/p- For completeness we include Qr, which is important in the 
early times. Taking time derivative of Eq.([2]) and making use of Eqs.([6]), ([7j), and ([8]) yield the 
following general form of the Statefinder pair 
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which hold actually for a generic dark energy model. 




Figure 2: The left panel: r{z) ~ 1 for z < 10 in the three cases of YMC. The right panel: s{z) ~ 
for z < 5. Thus (r, s) ~ (1, 0) with an error ~ 0.02 in YMC, quite close to that of LCDM. The initial 
YMC fractional energy il^j = 10~^ is taken. 

Let us apply the pair to some simple models. For any LCDM with a non-zero cosmological 
constant A, one simply has (r, s) = (1,0), and q{t) = (—1, |). For the Steady State Universe 
(SSU) model |29j| with a{t) = e^°^ one also has (r, s) = (1,0), but it has a fixed q(t) = — 1. 
For the standard cold dark matter model (SCDM), our calculation gives (r, s) = (1, |), which 
is different from the value (r, s) = (1,1) often quoted in literatures t22j. The details of our 
derivation is given in Appendix. 

In the following we present (r, s) in the YMC model for three cases: the non-coupling, the 
YMC decaying into matter, the YMC decaying into both matter and radiation. The following 
plots are based on Eqs.([TO]), (fTTD and (fT2|). 

First, we discuss the evolutions of (r, s). Fig|2] gives the behaviors of r{z) and s{z) in the 
recent past time with z < 10 in all the three cases of YMC model. It is clear that r{z) ~ 1 for 
z < 10 and s{z) ~ for z < 5. Thus (r, s) ~ (1, 0) with an error ~ 0.02 in YMC, quite close to 
that of LCDM. The two coupling cases, i.e., YMC decaying into matter, and into both matter 
and radiation, yield almost overlapping results since the YMC-radiation coupling T'/Hq ~ 10^"^ 
is very small. Furthermore, we would like to analyze the behavior in the radiation stage. In 
the limit of large z, one has ily ^lyi ~ (10^^°, 10~^) and 1 — ^lyi ^ 1, Eqs.l fTTl) and 

(fT2|) give an asymptotic behavior 
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independent of the initial condition. The analysis for the 1-loop YMC model in Ref. [T7] did 
not include Qr, so the extrapolation to the radiation stage could not have been made. 

To show the difference of YMC from other dark energy models. Fig 0] demonstrates r{z) 
and s{z) from YMC, LCDM, quiessence, and kinessence, respectively [22j. On the left panel, 
while YMC gives r{z) = 1 ±0.02 around 2; ~ 0, quiessence and kinessence give their respective 




Figure 3: The details of r(z) on the left and s{z) on the right around z = (0, 10) in YCM compared 
with quiessence (Q) and kinessence (K) models [22] . 

r{z) decreasing to 0.5 and 0.4. On the right panel, for 2; < 4 in all three cases of YMC model, 
s{z) approaches to 0, quite close to LCDM model. On the other hand, the quiessence has a 
constant s = 0.2, and the kinessence has s{z) > 0.3 around 2 ~ 0. The differences of the YMC 
from quiessence and kinessence are very large. However, like the other dark energy models, 
s{z) in YMC also has a divergence at 2; ~ 13, which does not show up in Fig 131 

In the following we discuss the r — s plan and the r — q plan in the three cases of YMC 
model. 

Figlll(a) gives the recent {z < 0.8) (r, s) trajectory in the non-coupling case of YMC. Two 
trajectories are presented for two initial values t/i = 1 and t/i = 15, respectively. The round dot 
at (1.00016,4.76x10-^) is the current value for = 1, and that at (0.997,9.07x10"^) is for = 
15, respectively. The star at (r, s) = (1, 0) is the fixed point of LCDM model, which currently 
is still robust against the observations. Our model predicts an asymptote at t ^ 00 very close 
to that of LCDM model. The trajectories of i/i = 15 and of yi = 1 approach to it by different 
paths. With the expansion of the universe, all trajectories for different initial conditions will 
reach to the fixed point (r, s) = (1, 0) ultimately. In comparison, the quiessence and kinessence 
models give the current values r < 0.5 and s ~ 0.5, other models have very scattered typical 
values, such as the quietessence model ^24j with (r, s) = (0.4,0.3), the Chaplygin gas model 
[25] with (r, s) = (1.95,-0.3), the agegraphic model [26] with (r, s) = (—0.2,0.5), which are 
far away from that of LCDM model [22] . The holographic model [27] without interaction has 
(r, s) ~ (0.94, 0.01), but an interaction 6^ = 0.1 gives (r, s) ~ (0.75, 0.09), deviating away from 
LCDM model again. 

Figlll(b) gives the recent [z < 1.1) (r, g) trajectory in the non-coupling case of YMC. For 
Hi = 1, r decreases to unity monotonically, whereas, for yi = 15, r decreases to a minimum 
value then increases back to unity. The two trajectories eventually approach to the fixed point 
at (1, —1) of the SSU. Note that, the pair (r, q) in YMC does not pass through the fixed point 




Figure 4: The recent trajectories (r, s) and (r, q) for the non-coupling case in (a) and (b) and for the YMC 
decaying into matter in (c) and (d). The round dots are the current values for the yi — I and for yi — 15, 
respectively. The arrows along the curves denote the evolution direction. The square dot at (5,1) is the fixed 
point of SCDM, and the triangle dot at (—1, 1) is that of the Steady State Universe. 

(l,i) of the SCDM pi]. 

Fig|4](c) shows the (r, s) trajectory in the case of YMC decaying into matter, whose overall 
profile looks similar to, but its detail is different from, that of the non-coupling case in FigJl](a). 
The current values of the Statefinder pair (r, s) = (0.981, 5.27 x 10^^) for yi = 1, and (r, s) = 
(0.979,6.21 X 10^^) for yi = 15, respectively. As t — > 00, (r, s) will also approach the fixed 
point (1, 0) of the LCDM model. 

Fig|l](d) shows the (r, q) trajectory in the case of YMC decaying into matter, which is 
similar to that of non-coupling case in Figj4](b), but the case for y^ = 1 has r < 1 due to 
the coupling. As t ^ 00, (r, g) will also approach the fixed point (1, —1) of the SSU models. 
Moreover, in this case, the trajectories of (r, q) do not pass through the fixed point(l, |) of the 
SCDM either. 

The case of YMC decaying into both matter and radiation. Their r{z) and s{z) are plotted 
in Figj2l Since the YMC-radiation coupling T'/Hq is very small, its modifications to the 
parameters r, s, and q are < 10~^. Therefore, the trajectories (r, s) and (r, q) in this case are 
almost overlapped with FigjH (c) and (d), respectively. To save room, we do not plot the (r, s) 
and (r, q) in this case. 

4 Summary 

The Statefinder pair (r, s) is examined for the 2-loop quantum effective YMC dark energy 
model. Three cases are presented: the non-coupling YMC, the YMC decaying into matter. 



and the YMC decaying into both matter and radiation. It is found that in all the three 
cases the pair (r, s) is very close to the fixed point (1, 0) of LCDM model for z < 5, and 
the deviations are tiny 6r ~ 10~^ and 6s ~ 10~^. Among the three cases, (r, s) in the non- 
coupling case differs only by ~ 1% from those in the two coupling models, while the two 
coupling models are almost the same as each other in all aspects, since the decay rate of 
YMC into radiation is very small in the model. In regards to the diagnosis of Statefinder 
pair, the YMC model is shown to differ drastically from other dark energy models, such as 
quiessence, kinessence, quintessence, Chaplygin gas, interacting holographic, and agegraphic, 
etc. If further cosmological observations continue to support LCDM model, they are unlikely 
to rule out the YMC dark energy model by using only the pair (r, s). 

ACKNOWLEDGMENT: Y.Zhang's research work was supported by the CNSF No.l0773009, 
SRFDP, and CAS. 

Appendix 

In this appendix we derive the Statefinder (r, s) for SCDM with A; = 0, i.e., the Einstein- 
de Sitter model containing only matter. For the time being, however, we let the radiation 
density pr be non-vanishing and will set pr = in the final step of calculation. The Friedman 
Equations are still given in Eqs.([l]) and ([2]), while the total energy density p = Pm + Pr and 
the total pressure p = Pr- There is no coupling between matter and radiation, so the energy is 
conserved for each component: 

pm + 3Hpm = 0; 

Pr + 3H{pr+Pr) =0, (14) 

which ensure that the total energy satisfies: 

p + 3H{p + p) = 0. (15) 

Taking time derivative of Eq.(l2]) leads to 
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Then with the help of Eqs.([l]), ([2]) and ( |T5l) . one obtains, 
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Applying Eq. (fT^ and taking p = Pr = \pr into account, one obtains r for SCDM model. 



r = 1 + 2nr. 
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By Eq. ([2]), the deceleration parameter q for SCDM model is given by 



" ^ = 37^ = TT. (20) 



Using the definition in Eq.Q, one has 

r - 1 _ 2^]^ _ 4 
3(g- 1/2) ~ p; ~ 3' 

which is diff'erent from the value s = 1 that has been often quoted in literature [22j. Note 
that, since VL^. is cancelled, this result of s is independent of the value of ^2^. Setting fi^ = in 
Eq.([T5|) yields the Statefinder pair of SCDM 

(r,s) = (l,^). (21) 
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